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Abstract. In this paper we discuss the notion of the diffusive capacity for 
the generalized Forchheimer flow of fluid through porous media. The diffusive 
capacity is an integral characteristic of the flow motivated by the engineering 
notion of the productivity index (PI), [81I1T10. The PI characterizes the well 
capacity with respect to drainage area of the well and in general is time depen¬ 
dent. We study its time dynamics for two types of fluids: slightly compressible 
and strongly compressible fluid (ideal gas). In case of the slightly compressible 
fluid the PI stabilizes in time to the specific value, determined by the so-called 
pseudo steady state solution, 013111 ]. Here we generalize our results from [3] 
on long term dynamics of the PI in case of arbitrary order of the nonlinearity 
of the flow. 

In this paper we study the mathematical model of the PI for compress¬ 
ible gas flow for the first time. In contrast to slightly compressible fluid this 
functional mathematically speaking is not time-invariant. At the same time it 
stays ’’almost” constant for a long period of time, but then it rapidly blows up 
as time approaches the certain critical value. This value depends on the initial 
data (initial reserves) of the reservoir. The “greater” are the initial reserves, 
the larger is this critical value. We present numerical and analytical results for 
the time asymptotic of the PI and its stability with respect to the initial data. 
Using comparison theorems for porous media equation from m we obtain 
estimates between the Pi’s for the original gas flow and auxiliary flow with 
a distributed source. The latter one generates the time independent PI, and 
can be calculated using formula similar to one in case of slightly compressible 
fluid. 
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1. Historical Remarks and Review of the Results 

The classical equation describing the fluid flow in porous media is the Darcy’s 
law, stating the linear relation between the pressure gradient Vp and velocity u. 
Darcy himself observed in |3] that the area of applicability of linear relation is 
very limited. When, for instance, the fluid has high velocity or in the presence of 
fractures in the media, the nonlinear models are necessary to capture the properties 
of the flow. 

One of the widely-used nonlinear models is the Forchheimer equation in the form 
g{\u\)u = —Vp, where g{s) is a polynomial (see [31 [T3]). Originally Forchheimer 
in his work [ 12 ] proposed three particular equations to match the experimental 
data: two term, three term and power laws, with g{s) being up to second order 
degree polynomial. To embrace the recent findings on the nonlinearity of the fluid 
flow (see mm) and simplify the mathematical handling, it is convenient to con¬ 
sider the generalization of classical Forchheimer equations to the case where g{s) is 
the generalized polynomial with non-negative coefficients and, possibly, non-integer 
powers (see m)- We call this family of equations the g-Forchheimer equations. The 
p-Forchheimer equation combined with the equation of state and the conservation 
of mass results in a single degenerate parabolic equation for the pressure p(a;,t), 
only. 

In this paper we discuss two types of fluids: slightly compressible fluid, char¬ 
acterized by the equation of state p(p) ~ exp( 7 p) with very small compressibility 
constant (7 ~ 10 “®), and strongly compressible fluid, ideal gas, characterized by the 
equation of state pip) ~ p, [18]. Note that in our model we assume that the porosity 
of the porous media does not depend on the pressure. In case of slightly compress¬ 
ible fluid we studied different properties of p-Forchheimer equations [HlllISI. In 
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this paper we extend the results of our work j3] on asymptotic behavior of the pres¬ 
sure function to the case when the degree of the ^-polynomial is arbitrary. In case 
of ideal gas we discuss both numerical and analytical results for the time dynamics 
of the solution of the corresponding parabolic equation. 

Keeping in mind that the applications of our findings are in geophysics and in 
particular in reservoir engineering, we restrict out studies to fluid flow in a reservior 
with bounded domain U. The reservoir is bounded by the exterior impermeable 
boundary Tg and the interior well-boundary T^. is subject to various boundary 
conditions depending on the flow regime. We introduce the capacity type functional 
to study the asymptotic behavior of the fluid flow in the reservoir with respect to 
time t. This functional is motivated by the notion of the well Productivity Index 
(PI) and is often used by reservoir engineers to measure well capacity, see [HI HI] . 
Productivity index is the total amount of fluid per unit pressure drawdown (the 
difference between reservoir and well pressures) that can be extracted by the well 
from a reservoir (see [HI]). It is defined as Jg{Ti){t) = Q{t)/PDD{t), where Q{t) 
is the total flux through the boundary T^, and the pressure drawdown PDD{t) is 
equal to the difference between averages of the pressure in the domain and on the 
boundary P^. 

In general the diffusive capacity is time dependent (see [HI mi) for both slightly 
compressible and compressible fluid flows. However its time dynamics differs greatly 
in these two cases. 

For slightly compressible fluids and specific regimes of production the PI stabi¬ 
lizes in time to a constant value. This value can be determined using the solution of 
a particular boundary value problem. Namely, the time-invariant diffusive capacity 
is associated to the pressure distribution ps{x,t) = - + W{x). Here the initial 

pressure distribution W{x) is the solution of the specific steady state BVP, Qs is 
a constant flux through the well-boundary, and V is the volume of the reservoir 
domain. Such pressure is called the pseudo-steady state (PSS) pressure and satis¬ 
fies the split boundary condition on the well —^t-\- p{x) for some known function 
(p{x). 

For arbitrary initial data one can not expect the PSS pressure distribution and 
constant diffusive capacity. However, as it appears from the engineering practice 
(see [HI mi) the constant PSS PI serves as an attractor for the transient PI. We 
proved this fact in mm under a number of conditions. In this paper we will 
generalize our previous results. We consider two types of boundary value problems. 

In the IBVP-I we impose total flux Q{t) on the well-boundary and consider the 
trace of the pressure function to be split as 7 (t) -|- tp{x,t), with /p tpds = 0. In 
this case 7 (t) can be considered as the average of the trace function and ip as the 
deviation of the actual trace from its average. Note that we impose conditions only 
on the function i/i, while 7 is unknown. For t —>■ 00 the boundary data {ip{x, <), Q{t)} 
are assumed to be localized in the neighborhood of the time independent {p{x), Qs}- 

In the IBVP-H no assumptions on the well-boundary flux are made. Instead, 
we specify the Dirichlet condition j{t) -I- tpixjt) on the well boundary. In this case 
both 7 and ip are known functions. For t —>■ 00 the boundary data {ip{x,t),j{t)} 
are assumed to be localized in the neighborhood of functions {ip{x),—^f-t}. 

The main result for both IBVP-I and H is that the time dependent diffusive 
capacity asymptotically stabilizes in the neighborhood of the diffusive capacity for 
the PSS regime associated with the pair (p{x) and Qs (see Theorems 4.11 and 5.4). 
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The corresponding value for the steady state diffusive capacity can be calculated 
just by solving an auxiliary Dirichlet time independent BVP. 

In [3] we obtained this result under the degree condition, stating that the degree 
of ( 7 -polynomial deg(( 7 ) < 77 ;^, where n is the dimension of space. Mathematically 
this condition arises from the theory of Sobolev spaces and insures the continuous 
embedding C L‘^{U). While for n = 2 this constraint holds for any 

degree deg{g), for n = 3 it holds only for deg(( 7 ) < 4. In this paper we study the 
asymptotic behavior of the diffusive capacity without any constraint on deg{g). In 
this case the classical Poincare-Sobolev inequality doesn’t hold, and we use weighted 
inequality for the mixed term a ^ 2. We obtain estimates for the 

bounds of the L“-norm for both the difference between transient and PSS solutions, 
and the difference between transient and PSS solution time derivatives. 

In Sec. we discuss the concept of the diffusive capacity for an ideal gas flow 
and some results on its time dynamics. In this case the productivity index is 
defined according to [317] as Jcit) = Q{t)/PDD{t), where PDD{t) is the difference 
between the average of in the domain U and on the boundary P^. In contrast to 
the case of slightly compressible flow, we show numerically that until certain critical 
time Tcrit the transient PI remains almost constant and then as time approaches 
Tcrit it blows up. This result obtained on actual field data corresponds with the 
engineering observations. Time Tcrit depends on the initial reserves/initial data: 
the “greater” the initial reserves are the greater Tcrit is. Similar to our approach in 
case of slightly compressible fluid, we use the auxiliary pressure poix, t), resulting in 
time independent PI to investigate the behavior of time dependent PI. The Po(x, t) is 
the solution of the equation with positive function on the RHS. This function can be 
considered as fluid injection inside the reservoir. When gas reserves are considerably 
larger than pressure drawdown on the well, t his source term is negligible, and the 
Pi’s for po{x,t) is almost identical with the general time dependent PI. For linear 
Darcy case we obtain some analytical comparison theorems between the pressures 
for time t < Tcrit- Under some constraints on the boundary data and smoothness 
of the pressure function we obtain a stability result for the PI with respect to the 
initial data. 

The paper is organized as follows. In Sec. we discuss the various aspects of 
nonlinear Forchheimer equations and associated parabolic equation. We give the 
definition of the diffusive capacity on the solution of this parabolic equation. The 
corresponding IBVP-I and II for the total flux and Dirichlet boundary conditions are 
introduced in Sec.|^ Sec.j^is devoted to the IBVP-I: we state the conditions (4.7) 
- (4.10) imposed on the boundary data to obtain the main result on the asymptotic 


convergence of the transient diffusive capacity (see Theorem 4.11). In Lemma 4.1 


we specify the estimates on the pressure, its gradient and time derivative that are 
necessary to prove this result. Sec . [5] i s devoted to the similar results for IBVP- 
II: we state the conditions (5.4) - (5.8) imposed on the boundary data to obtain 
the asymptotic convergence of the transient diffusive capacity (see Theorem 5.4). 
Finally, in Sec. |^we discuss the concept of PI for an ideal gas flow. 


sec:problem-state 


2. Problem statement and Preliminary Properties 

Consider a fluid in a porous medium occupying a bounded domain U C K". Let 
X G M", n > 2, be a spatial variable and < S K be a time variable. Let u{x, t) G K" 
be the fluid velocity and p{x, t) S K be the pressure. 
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eq:g-forcli 


non-lin_darcy 


Kg 


eq:continuity (2.6) 


We consider a generalized Forchheimer equation 
( 2 . 1 ) g{\u\)u = -Vp, 

where g : M"*" —>■ K’*'. In particular we consider function g to be the generalized 
polynomial with non-negative coefficients. Namely 

k 


gdef (2.2) 


g{s) = aos“° -f ais“^ H--I- Ofcs“'‘ = oo + ^ ajS° 




with A: > 0, the real exponents satisfy ao = 0 < ai < a 2 <■■■< cxk, and the 
coefficients ao,ai ,... ,ak > 0. The largest exponent is the degree of g and is 
denoted by deg{g). 

One can notice that the equation (O) with g{s) defined as in ( |2.2[ ) includes the 
linear Darcy’s equation and all classical forms of Forchheimer equations [l2] (for 
details see our previous works mm- 

In case of slightly compressible fluid we consider the case when the Degree Con¬ 
dition in (2.14), [3], is not satisfied, namely 


ndeg-cond (2.3) 


Ufc 


< 


Uk = deg(ci) > 


afe-|- 12 -|-n n — 2 

In this case there is no continuous embedding C Lp'iU) and correspond¬ 

ing Poincare inequality doesn’t hold. Clearly, if n = 3 condition (2.31 will hold for 
the nonlinearities ak > 4. 

From (2.1) one can obtain the non-linear Darcy equation explicitly solved for 

u = —KQVpDVp, 


the velocity u: 

(2.4) 
where 

(2.5) K{0 = 


1 


^ > 0, G(s) = sg{s), s >0. 


9{G-^{oy 

Along with which is considered as a momentum equation, the dynamical 

system is subject to the continuity equation 


| + v.w = o, 


where p{x, t) G M"*" is density of the fluid. The fluid is considered to be slightly 
compressible, subject to the equation of state 


eq:state (2.7) 


^ = 7 P or p(p) = poe^^P 


where 7 ^ 10 ® is the compressiblity constant, see [Ill- 

Combining (2.1), (2.6) and (2.7) we obtain the degenerate parabolic equation 
for pressure (see [ 1 ] for details): 


eq:p-l (2.8) 


7 ^=V-(iF(|Vp|)Vp). 


By scaling the time variable we can assume further on that 7 = 1 . 

In case of the flow of ideal gas the equation (2.1) for n = 2 will take the form 
(see, for example, [13] ') 

au +/3pu\u\ =—Vp, /3 = 
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eq:gas-intro (2.9) 


As in case of slightly compressible fluid, the above equation can be solved for u: 

2 


u = — 


a + \/a^ + 4/3p|Vp| 


Vp. 


Combining the equation above with the equation of state p{p) = p and continuity 
equation (2.6 1 we obtain (for the details see Sec. 


dt 


2p 


a + 4/3p|Vp| 


Vp . 


def:av-pressure 


We will study equations (2.8) and (2.9) in the open domain U C K" with the 


boundary dU = F = Fe U F^. The F^ is considered to be external impermeable 
boundary of the reservoir with u ■ N = 0, where N is the outward normal vector 
to Fg. The Fi is the internal boundary of the well with the total flux or Dirichlet 
boundary conditions. 

To study the long term dynamics of the i^-Forchheimer flow in the domain U 
we introduce special capacity type functional. We call it the Diffusive Capacity/ 
Productivity Index (PI) (see [HI HIl HH] ) and define as 


def:pi (2.10) 


JgW — 


Q{t) 

PDD{t) ■ 


where Q{t) is the total flux through the well-boundary F^ and PDD{t) is called the 
pressure drawdown in the domain U. 

In case of slightly compressible fluid these quantities are deifned as 


( 2 . 11 ) 


Q{t) = 


u ■ N ds, 


(2.12) PDD{t) = p^{t) - py^(t) = j^p{x,t)dx - p{x,t)ds. 

are the pressure averages in the domain and on the well-boundary correspondingly 
For gas filtration in porous media these quantities are deifned as, see 00: 


pu • N ds, 


Q{t) = ^ 


We end this section with the recollection of some properties of Ar(^) that will be 
used in this study (see [I||3]). Function K{^), ^ G [0, oo) is decreasing and satisfies 


k 

-prime 




k-ineq 


k-xi2 


(2.13) 

(2.14) 

(2.15) 


0>if'(?)>-a^, 

Co(l + ?)■“ < < Ci(l + 0”“ < Cl, 

C'2^2-a _ ^ 
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ineq:monoton 


def:c-phi 


h-def 


h-k-ineq 


h-xi-2-a 


sec:ibvp-dc-shift 


eq:p 


be:flux 


be: O-neumaoi 


initialcond 


trace-p 


Moreover the function K{\y\)y, y G K", is monotone, i.e. for any two functions 
ui,U 2 € one has 

(2.16) / (X(|Vwi|)Vui-/ s:(|Vm 2 |)Vu 2 ) • V(wi 
Ju 

> |V(mi - U 2 )\^~°‘ dx 

where 

(2.17) (7$ = C'$(ui,M2) = Ci{l + max(||Vui||i2-a(y), |jV w2||l2-“(c/)))‘ 
Following mm we define the functional 

(2.18) iJ(C = [ K{yfs) ds for ^ > 0 . 

Jo 

Function H{^) can be compared with K{^) (see [T]) as follows: 

(2.19) K{Oe < HiO < ‘^K{i)e. 


Combining (2.191 with (2.15) there exist constants Ci and C 2 such that 
( 2 . 20 ) - 1 < H{C) < ( 72 ^^”“. 

3 . Two TYPES OF BOUNDARY CONDITIONS AND PSEUDO-StEADY StATE 

SOLUTION 

In Sections and we consider the case of slightly compressible fluid. 

Let U C K" be an open set with the (7^ boundary dU = F = Fg U Fj. The Fg is 
considered to be external impermeable boundary of the reservoir and the F^ is the 
internal boundary of the reservoir (well surface). With respect to boundary F^ two 
different problems will be considered: IBVP-I with imposed total flux (generalized 
Neumann condition) and IBVP-II with imposed Dirichlet boundary condition. 

I. IBVP-I for the total flux condition. The function p{x,t) is a solution of the 
IBVP-I if it satisfies 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


Qt) 

^ = V • (77(1 Vp|)Vp) in D = Ux{0, 00 ), 

— J K{\Vp\)Vp ■ N ds = Q{t) on rix( 0 ,c»), 


dp 


= 0 on Fe X (0, 00), 


p{x,0)=po{x) in U. 


Since the IBVP-I ( |3.1[ )-(3.4) lacks uniqueness, following [3] we restrict the boundary 
data. In particular, let the trace of p{x,t) on well-boundary F^ be of the form 


(3.5) 
with 

(3.6) 


P{x,t)\ri =7(t)-b'0(ai,t), 


J ifix, t) ds = 0. 


Note, that function 7 ( 7 ) is not specified and will be determined by the flux Q{t). 
Such restriction of the boundary trace ensures the uniqueness of the solution of 
IBVP (3.1)-(3.4) (see Remark 3.2 in [3]). 
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be:dir 


deviation-dir 


pss-profile 


eq:w 
be :w 


be:w-neum 


def:pi-pss 


Transient boundary data tjj{x,t) and Q{t) will be compared with the time- 
independent boundary data (p{x) and Qs- As in [3] we consider functions 
and <i)(a:) to be defined in the domain U with the traces on being t) and ip{x) 
correspondingly. The results are obtained under the constraints on the parameters 
of difference of boundary data 

(3.7) Agit) = Q{t) - Qs] Aq,{x,t) ='i>{x,t) - ^{x). 


II. IBVP-II for Dirichlet boundary condition. The function p(a:, f) is a solution of 


the IBVP-II if it satishes (3.11, (3.3), ( |3.4[ ) with the condition on the well boundary 

(3.8) p{x,t)\Ti = ii{x,t)+^{t) on rjx(0, oo), 


for known function ip{x,t). 

Transient boundary data %f{x^t) and ^{t) will be compared with the boundary 
data p{x) and —At, A = Qs/\U\, where Qs is constant flux on T^. We again 
consider functions 4'(a;,t) and $(a;) to be the Wf{U) extensions of if{x,t) and p{x) 
on the domain U. The results are obtained under the constraints on the differences 


(3.9) A\i,{x,t) ='i{x,t) — ^{x)] A.y {t) = ^{t) + At. 

For the existence and regularity theory of degenerate parabolic equations, see 
e.g. [Hiniin] and references therein. 


3.1. PSS solution. For the time-independent boundary data (p{x), Qs we define 
the pseudo steady state (PSS) solution of IBVP for equation 


(3.10) 


Ps{x,t) = -At + W{x), 


\u\ 


for a given function (p(x) and a constant total flux Qs on the boundary F^. Function 
W{x) is called the basic profile corresponding to the flux Qs and is defined as a 
solution of BVP 


(3.11) 

(3.12) 

(3.13) 


-A = V - {K{\VW\)VW), 
Wix)\ri = (fix), 
dW 


dN 


= 0 . 


The existence and uniqueness of a solution for the given pair (p{x), Qs is proved in 

0 . 

It is not difficult to see that in case of the PSS solution when p{x,t) = ps{x,t) 
the functional Jg(t) in (2.10) is time independent and 


(3.14) 


Jg{t) = Jg,PSS = 


Qs 


\u\ fu ^(^) ir^i /p; 


= const. 


We assume that <p{x) and Qs are such that the denominator of (3.14) is not equal 
zero. Obviously PSS solution is the only solution for given (p{x), Qs which obey 
both types of boundary conditions. Therefore PSS will serve us a reference solution 
for comparison for both IBVP-I and II. 
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thrm:main-gen 


lem:pi-diff 


cond:j-js 


J-JPSS 


3.2. Time convergence o f Jq {t) Jg,pss - general case. Let p{x,t) be the 
transient solution of IBVP-I (3.1 )-(3.4) with boundary data 4'(a;, t), Q(t) or IBVP-II 
(3.1), (3.3), (3.4), (3.8) with boundary data Let Jg{t) be the corresponding 

diffusive capacity defined as in (2.10). Letps(a;,t) be the PSS solution of this IBVP 
with boundary data ^{x),Qs and the corresponding constant diffusive capacity 
Jg,PSS defined as in (3.14). 

For slightly compressible fluid flow subjected to Froschheimer equation the aim 
of this study can be formulated as follows 

Theorem 3.1. There is a wide class of transient boundary data such that for both 
problem IBVP-I and IBVP-II 

~ Jg, PSS 0 as t —>■ 00 . 

The proof uses the following Lemma which reduces the question of the conver¬ 
gence of J{t) to the convergence of the gradient of the solution in the appropriate 
norm and of the total fluxes to a corresponding constant values. 

Lemma 3.2. As t > oo the difference between transient and PSS PPs 

(3.15) Jg{t) - Jg^pss ^ 0 if \\^{p - Ps)\\l^-^{u) ^ 0 and Aqf) ^ 0. 

Proof. From (2.10) and ( |3.14 ) it follows 

Q(t) Qs Q{t){PsU-P stJ - QsiPu-PrJ 


Jg(t) - Jg.PSS = = 


Pu-Pr, PsU-PsT, (Pu -PrJiPsU -PsT.) 

Q{t)[{PsU ~ PsTj) ~ (Pu ~PPj)] + ^q{^){Pu ~PrJ 
{Pu -PtMPsU -PsT,) 


(3.16) 


Q{t) ■ Ap{t) 




J, 


g,PSS 


(Pu - PrJiPsU - PsF.) ' Qs 

Here Ap{t) = (p^p — p^p.) — (pjj —Pr ) is the difference of pressure drawdowns for 
transient and PSS regimes. 

The pressure drawdown for the PSS solution is time independent, and is not 
equal zero. Hence Jg{t) — Jg,pss —>■ 0 if and only if both Ap{t) — >• 0 and Agf) — >• 0 
as t —>■ oo. 

The difference in pressure drawdowns Ap{t) can be written as 


Ap(t) - 1^, 
1 


J {p - Ps) ds - J {p - Ps) dx 


|rz 


= THT / [ P - Ps - Jm ip - Ps) dx ] ds 


\U\ 


1 

W\ 


P-Ps- 


\u\ 


jjyP- Ps) ^ |I^ ^ ^ 


where 


z{x, t) = p - Ps - j ip-Ps) dx. 

Applying Trace theorem and Poincare inequality, as fjjzdx = 0, one has 
|Ap(t)|2-“ < C / m^-‘^ds<C [ |Vzp-“cix = C' / |V(p - p,)p-“ dx. 

OFi JU Ju 

The condition (3.15) follows. 


□ 
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sec:flux 


def:fl 


def:f2 


The results for the solution of IBVP-I will be obtained under the conditions on 
the deviation of the boundary dat a (|3.7[ ) in Sec. and for the solution of IBVP-II 
in terms of deviations (|3.9|) in Sec. [Sj 


4. IBVP-I: Asymptotic convergence of Diffusive Capacity 


In this section we prove the Theorem 3.1 on asymptotic convergence of transient 
diffusive capacity Jg{t) defined on the solution of IBVP-I (3.1)-(3.4) with boundary 
data 4'(a;, t), Q(t). Let ps{x,t) be the PSS solution of this IBVP with boundary 
data 4>(x),(5s and the corresponding constant diffusive capacity Jg^pss- 

4.1. Assumptions on boundary data. In order to formulate the assumptions 
in the Theorem |3.1| we first introduce the following notations. 

Suppose a > > 2 and let 2 — a < 6 = < a. We define 

(4.1) Fi(t) = AQ(T)dT+J^Aq,dx-J^(p(x,0)-ps(x,0))dx. 

(4.2) F2(t) = 1 + HVWh2-.(p) + lF,(t)l + f |V(A^)|dT+|AQ(t)|T^. 

Ju 


def:al 


def:limsupal 


def:a2 


def:a3 


(4.3) 

(4.4) 

(4.5) 


Ai(a,t) = IIVWII 


OL — a 


1 ^* 1 “ dT 


+ llv(A>,)ii“r 


-b 




it). 


lb/ dx 


Al (a) = lim sup Ai (a, t) . 

t—¥oo 

A2{t)= [ \\7W\^-^dx+ [ iV^'tpdx-b [ |V(A^)|2-“£ia; 
Ju Ju Ju 

\'^t\^dx 


\'^t\dx 


2-a 


\'^t?dx- 


2 —a 


+ |AQ(t)| + \Q'{t)\ + |AQ(t)|^ + \Q\tt + C\Fi{t)\ 

(4.6) A3{e,t)= [ \W'$Udx+i-^ [ \^tUdx+\Q'{t)f + \F[{t)f. 

Ju Ju 

The certain way the conditions A1-A5 will be formulated is dictated by the fact, 
that while the main result will be obtained under all the assumptions below, some 
intermidiate inequalities require less restrictive constraints. As a note for each as¬ 
sumption we state the resulting bounds on parameters Fi, F 2 , Ai, A 2 and A 3 . 


assumpl 


Assumptions 1 (Al): (Ai(a) < C) 
(4.7) limsup (|AQ(t)|-b ||4't||L6 

t—^OO ^ 


IV4-II 


L’’ 


|V$|! 


Lb ^ 


< c. 


Assumptions 2 (A2): (Ai(q)- b limsup(A 2 (t)-b P 2 (t)) < C) 


assump2 (4.8) 


Assumptions Al; 
limsup (iQ'(t)l -b |Ai(t)| 


IV4' 


t||L= ^ 


< c. 
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assumpS-beta 


assump4 


assumpS 


def :q 


def:qs 


eq:q 


be :q 


eq: qs 


be: qs 


intq=0 


lem:nppsO 


eond:sometoO 


Assumptions 3-/3 (A3-/3): for any /3 > 2 

/ OO 

dt < C. 

Assumptions 4 (A4): limsup43(t) = 0 

t—^OO 



[ Assumptions A2; 


] limsup(^(|V4/,p + |vIi,,n 

t—^OO ^ 

1 dx + Q'(t)P + ^Q{t)j 


= 0 . 


Assumptions A3 and A4; 
limsup ||V(A^)|| 2,2 = 0 . 


Along with the solution of IBVP-I p{x, t) and the PSS solution Ps{x, t) of IBVP-I 
with boundary data $(a;), Qs, we will use the following shifts of the solutions: 

(4.12) q{x,t) =p{x,t) - J^p{x,t) dx - (^'i’{x,t) - J^'i/{x,t) dx'^ ; 

(4.13) qs{x) = ps{x,t) - J ps{x,t)dx - (^^{x) - J ^(a:)^^^ 

It follows that q satisfies the BVP 

(4.14) qt{x,t) =pt{x,t) + J^'^t{x,t) dx'^ , 

(4.15) q\r, ='y{t) - J^p{x,0)dx +J Q{T)dT +J^'S{x,t) dx. 


Similarly qs satisfies the BVP 
(4.16) 


dqs, .. 9ps 

^(i,() = ^(x,() + /l = 0, 


(4.17) 

Note that 

(4.18) 


Qs\r, = -|^ 

/ q{x,t)dx = 0 and / qs{x,t)dx = 0. 

Ju Ju 


According to Lemma 3.2 we need the convergence ^Q{t) 0 and ||V(p — 

Ps)||L2-a —>■ 0 as t —> OO. The affinity of Q{t) and Qs will be imposed as a condition 
on the boundary data (see Assumptions A4). The estimates necessary to prove the 
convergence of ||V(p—Ps)||L 2 -a are outlined in the following Lemma. The references 
to the corresponding results further in the paper are given. 


Lemma 4.1. Let q and qs be defined as in (4.121 and (4.13). Then 
||V(p - Ps)||i2-a 0 as t-> OO 


*/ 

(4.19) lim (||V(A^)|U 2 + |AQ(t)||Ai(t) + 1|) = 0 
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and 


ineq:zah-bound (4.20) 


1. / \q — Qsf dx < C / \q — qs\°‘dx < C, see Theorem 

Ju Ju 


ineq:np-bound (4.21) 2 . / |Vpp °‘dx<C, see Theorem \4- ?| 


ineq:pta-0 (4.22) 


ineq:nps-bound (4.23) 


3- / \pt + A\^dx= / \pt — Ps.t\^ dx ^ Q, see Theorem\4-10 
Ju Ju 

4. j |Vpsp““dx= j dx < C, see Lemma 4-1, 

Ju Ju 


Proof. According to the inequality (5.20) in proof of Theorem 5.6 in [3] we have: 
n-p-ps-1 (4.24) ||V(p-ps)||| 2 -a < (7$ / \pt + A\\q - qs\dx 

Ju 


+ |V(A^)|2 dx + Alit) + \F,{t)AQ{t)\ + ||V(A^)|U 2 ^ . 


Here (7$ = C<i,{p,W) is as in (2.171 and is bounded by virtue of (4.21) and ( 4.23[ ). 
It is then clear that the RHS of (4.24) converges to 0 under conditions (4.19)- 


(4.23). 


□ 


4.2. Bounds for the solutions. As before let q{x,t) be the solution of BVP 
(4.14), (4.15) and qs{x,t) be the solution of PSS BVP (4.16), (4.17). Let z = q — qs- 
We will prove that under certain condition |z|“ da; is bounded at time infinity. 

In order to obtain the suitable differential inequality we will use the following 
result (see Lemma 2.3, m) 

Lemma 4.2. Let U be an open bounded domain in K" and a > If u\gu = 0, 
there exists constant C such that 

f K{\\'u\)\\'u\^\u\°‘-‘^dx>C f |Vup-“|u|“-2-Cdx 
Ju Ju 


ineq:w-poiiicare-simple-0 (4.25) 


> c 


\u\°‘ dx 


-c, 


where 70 = a/{a — a). 

We first prove the following lemma to estimate the boundary data. 


lemma:bdd-f or-q-qs 


Lemma 4.3. Let Fi{t) be as in (A.l), Aq = Q{t) — Qs and A.y = j(t) + At. Then 
/ora>^ 

(4.26) |A^(t)+Fi(t)r-i|AQ(t)| 


< Cl 


i 


'' dx 


-£2 / + 


with 7 o = a/{a — a). 
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Proof. We have zlr^ = {q — 9 s)ri = Then by Trace theorem 

(4.27) ((A^ + 7^i)-|r,|)“-“= ^ <cj^ 121“-“ ds 


= C < C / | 2 |“-“da; + C' / 

JTi Ju Ju ‘- 

|Vz||z|^l^ “ dx 




dx 


= C / | 2 |“-“da; + C' 

JU Ju ^ 


= C [ \z\°‘-‘^dx + C [ |V2|2-“|z|“-2da;. 
Ju Ju 


Then 


(4.28) |A^ + Fi|“-i <C\ j 121 


-1 , Q -1 

dx] “ + C f / |V 2 | 2 -“UI“-^cix 

-1 , . Q-1 

2 —a I ^ la —2 


<C[ I \z\‘^dxj |V 2 r-“| 2 |“-"dx 


Applying Young’s inequality with e we get 


|A^+Fi|“ ^|AQ(t)|<ei{ I \z\°dxj + J |V 2 | 


2 -“| 2 |“- 2 da; 


+ (Cei + Ce2)^o " ’ 


which proves the result. 


□ 


lem:ddt-q-qs-l Lemma 4.4. Suppose a > ffff- > 2 and let 70 = —Then for z = q — Qs and 


Ai{a,f) is as in (4.3) 


est :ddt-zbEir 


(4.29) 


^ J | 2 |“da;<-C'f / \z\°‘dx") + C{1 + Ai{a,t))- 


Proof. Subtracting equations (4.14) and (4.16) from each other, multiplying on 
| 2 |“-^sign( 2 ) and integrating over U one has 


1 d 

Q. ill J jj 


4 - f |2|“dx =-(a- 1 ) [ ( 7 s:(|Vp|)Vp-A:(|VW|)VW)V2|2|“-2d2; 

J't Ju Ju 

+ [ {K(\Wp\)Wp-K{\WW\)WW)-N\z\'^-^ds+ 4 -rAQ{t) [ |2|“-isign(2)dx 

JTi |C| JU 

~ Iu^'^*~ W\Iu^' k|“-isign(2) dx. 
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ddt-z-0 


4-eq 


Ip-final 


Ipw 


From boundary conditions (4.15) and (4.17) it follows that 
and one has 

-w / \z\°‘dx = -{a-l) [ iK{\Vp\)Vp- K{\VW\)VW)V{p- dx 

a dt Jjj Jjj 

+ (a-l) [ {K{\\7p\)\/p- K{\\/W\)\/W)\7 {A^)\z\°‘-^ dx 

Ju 


(4.30) 


1 

M 


[AQ(t)+ /" dx - |t/|^'t] |z|“ ^sign{z) dx + [A^{t) + Fi{t)]°‘ 

Ju 

We split the first integral in RHS of (4.30) in four separate integrals and estimate 
them one-by-one. Namely, 

(4.31) - [ {K{\\7p\)Vp-Ki\VW\)\7W)\7ip-p,)\z\‘^-^ dx = Ip + I^ + Ip^ + I^p, 
Ju 

where (notice that Vps = VIF) 

Ip = - [ if(|Vp|)|Vpnz|“-2da;, I^ = - [ K{\VW\)\\/W\^\z\^-^dx<0, 

Ju Ju 

Ip^= [ iF(|Vp|)(Vp-VW)|zr-2dx, I^p= [ iF(|VlF|)(VW-Vp)|z|“-2dx. 

Ju Ju 


Since Vz = V{p — W) — V(A,i,) and similar to (4.5) in [13] we 

have 


2—01 ^la —2 


Ip<-C [ (|Vpp-“ - l)|z|“-2dx <-c [ |Vz| 

Ju Ju 

+ C [ (|V(Ag.)|2-“ + + l)|zr-2 dx 

Ju 


dx 


<—CJ |Vzf “|z|“ ^dx + ei\ I |z|“cix 


(4.32) + Ce, [l + ||VW||“7-(“^) + ||V(A^)||“r(“^^_ . 

The integral Ipw can be estimated similar to (4.6)-(4.9) in [13|: 

(4.33) 

Ipw< f Ar(|Vp|)|Vp||VlF||z|“-2dx < [ |Vp|^-“|VlF||zr-2rfx 

Ju Ju 

< [ |Vz|^-“|VlF||z|“-2dx+ [ (|VlFp-“ + |V(A,t)|i-“|VlF|)|z|“-2(;a.. 

Ju Ju 

The first integral in ( 4.33| ) can be estimated similar to (4.8) 

[ |Vz|i-“|VlF||z|“-^dx 
Ju 


(4.34) 


<ey^|Vz|^-“|zr-^dx + e[ jJA^dx) +Ce\\VW\\lI^jjy 


The second integral in (4.33) is estimated similar to (4.7) 
(4.35) / |V1F| 




<e{ \zrdx\ +C,\\VW\\li;^^y 
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int-l-a 


Ipw-final 


Iwp 


Iwp-final 


Similar the third integral in (4.33) can be estimated as 
(4.36) f < e 


dx 


+ Cel^J \V{A^)\^\VW\- dx 
Since —I" = 1; from Holder and Young’s inequalities it follows 


(i-°) , 




(4.37) l^\\7{A^)\^\\7W\^dxj ' < C (||VlY||“7-(“^) + || V(A,.)||“7(t;)) ■ 


Combining the estimates above in (4.33) we have the estimate for the integral Ip^ 

1 

^ 70 

izr dx 


Ipw < e / |V0f “|z|“ ^dx + e 


'u 


'u 


(4.38) + C + ||V(A^)||2E(“^)) . 

Finally, for the integral ly^p we have 

(4.39) 

Iwp< f KilVWDlVWWVpWzf-^ dx < [ \VW\^-^\Vp\\z\°‘-^ dx 

Ju Ju 

<C [ \VW\^-‘^\Vz\\z\°‘-'^dx + C [ |VfF|i-“|V(VF +A,i,)||zr-2da;. 
JU Ju 


Applying Holder’s and Young’s inequalities for the first integral in the RHS of (4.39) 
we get 


[ |VVF|i-“|Vz||zr-^dx 
Ju 


(°-2) , . (c-2)(l-a) 


= / |VtY|i-“ • iVzll^l (5^ • |z |——dx 


< 


|Vz| 


^-^'zr-^dx 


|VVF| 




1 - Q 

2- a 


<s f iVzp-^l^r-^dx + Ce / |VVF|2-“|z|“-2da: 

Ju Ju 

<ef l^zl^-^^lz^-^dx + ef f +C,||VfF||“E(y). 


The second integral in the RHS of (4.39) can be estimated similar to (4.36) 

f |VfF|i-“|V(fF +A>f)||z|“-2dx 

Ju 

<e(f Iz^dxY^ +C'(||VfF||“E(t,) + l|V(A^)||“r(“^)). 


Then in (4.39) the integral /u,p can be estimated as 
(4.40) 


Up — ^ I IV z 


2-a|z|a-2 


dx + 6 |zr dx^ + C (||VfF||“7-(“y) + II V(A 


Y—a \ 

'*'hlLi>(U )j ■ 
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Ist-int-rhs 


2 nd-iiit-rhs 


some-bdd-est 


Substiting (4.32), (4.381 and (4.40) in (4.31) we get the estimate for the first integral 
in the RHS of (|430|) 


- / {K{\\/p\)\/p-K[\\/W\)\/W)\/{p-ps)\z\°‘~‘^ dx 

Ju 


<-{C-e) / |Vzr-“|z|“-"da; + e / \z\‘^ dx 


(4.41) 


C 


l + ||VTR||“r“^) + ||V(A^)||“r(t,) 


The second integral in the RHS of (4.30) can be estimated as 

(4.42) 

I {K{\\/p\)\7p - R:(|VlT|)VlT)V(Ag.)|z|“-2 dx 

Ju 

<C f K{\Vp\)\Vp\\V{A^)\\z\°‘-^ dx + f K{\VW\)\VW\\ViA^)\\z\°‘-'^ dx. 
Ju Ju 

Then similar to the estimate for and ( |4.36| ) we get 

[ {K{\\/p\)\/p-K{\\/W\)\7W)\7{A^)\z\°‘-^dx <e [ \\7z\^-^\z\°‘-^ dx 
Ju Ju 

1 

(4.43) + e |zr dx^ + C {\\^WrU^u) + II V(A^)||“r(t,)) ■ 

By Holder’s and Young’s inequalities, we get 

(4.44) ^ f [AQ{t)+ [ '^tdx-\U\-^t]\z\°‘~^sigiiiz)dx 

\^\Ju Ju 


<s{ Iz^dx] TC 


<e / Iz^dx] TC 


u 

AQ(t) 


AQ{t)+ / 'i’tdx- \U\ 


dx 


'll# dx 


+ C I'fA^dx 






Finally combining (4.41), ( 4.43[ ), ( |4.44[ ) and Lemma 4.3 in (4.30) we get after 
choosing small enough e 


yet-another-formula (4.45) 


jU Izl^^dx^-C |Vzp-“|z|“-2dx + e ( / \z\‘^dx] ‘ + (7(1 + Hi(a,t)). 


Applying (4.25) to the first term in the RHS of (4.45) we get 


d 

dt 


\zrdx < -C 


z^dx] -1 +e 


Selecting sufficiently small e we obtain (4.29). 


| 2 :|“fia;) + C{1 + Ai{a,t)). 

□ 


lem:Al-gronw 


In order to prove the estimate on jjj |z|“ dx we recall the Lemma A.l from [13]. 

Lemma 4.5. Let (j) : [0,oo) —>■ [0, c») be a continuous, strictly increasing function. 
Suppose y(t) >0 is a continuous function on [0, oo) such that 

y' < -C(t)~^{y{t)) + f{t), t > 0, 
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th:zbar-toO 


est:zbar-bound-2 


est:zbar-bound-1 


thrm:bdd-h-gradp 


ip-limsup 


ip-bounded 


qt-qst-6 


ilp 


eq:ddt-qqs-equal 


qt-qst-7 


qt-qst-8 


where f{t)>0 for t > 0 is a continuous function. Then 


limsupi/(t) < (p 

t—¥oo 


I limsup 

\ t—>-oo 


m\ 

h{t)) 


The following result will follow from Lemma [4.4| 


Theorem 4.6. For a > naf{2 — a) and Ai{a) as in (4.4) 

(4.46) limsup j \z\°‘dx < C Ai[a)). 

t—^OO J (J 

Consequently, if Q(t) and 'i>{x,t) satisfy condition A1 then 

(4.47) limsup [ \z\°‘dx < C. 

t^OO Jlf 

Proof. The result follows from Lemma [4. 5 1 


□ 


4.3. Bounds for the gradient of solutions. We will obtain the bounds for the 
integral H{\\^p\) dx, where iL(|Vp|) is defined in (2.18). This result is necessary 


for our estimates of the difference in time derivative of fully transient and PSS 
solutions. 


Theorem 4.7. Let A 2 {t) be as in (4.5). For a > na/{2 — a) 

(4.48) limsup j i7(|Vp|) dx < C ( 1 + 7li(a) + limsupy42(t)) . 

t—¥00 J If \ J 

Consequently, if Q(t) and 4!{x,t) satisfy conditions A2 then 

(4.49) J |Vp|2-“dx < C H{\Wp\)dx+l^ <C. 
Proof. In formula (4.46) in [5] we have ior z = q — qs'. 




h [p] (t) + z'^ dx 

Ju 


< — C Ji [p] (t) + C(1 + A2(t)), 


(4.50) 
where 

(4.51) Ii[p]{t)= I H{\\7p\)dx- I K{\\7W\)\/W\/pdx + AQ{t)A^{t). 

Ju Ju 

Notice that 

(4.52) — f z^ dx = 2 f zzt dx > — f z^ dx — f q^ dx. 

dt Ju Ju Ju Ju 

Subtracting (4.52) from ( |4.50 ) we get 

(4.53) ^ Ii[p]{t) <-CIi[p]{t) + J z'^dx + C{l + A 2 {t)). 

The term dx can be estimated using (4.46) as a > 2. Then in (4.53) we have 


(4.54) 


^ di[p\it) < —CIi[p\{t) + C [1 + A2{t) + Ai{a)]. 
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Applying Gronwall’s inequality one has for t > 0 


qt-qst-9 


h[p]it)<e Ii[p]{0) + C f {I + A 2 {t) + Ai{a)) dr 

Jo 

(4.55) <C(l + Ai(a))+e-"i*/i[p](0) + C / A^ir) dr. 

Jo 

According to formula (2.35) in [T3] 

est:limsup-int-Ot (4.56) limsup / e“^^*“’^^/(r) dr < limsup/c“^/(t). 

- t —¥00 Jq t —>00 


estim-H-final 




ddt-q-qs-5 




dgamma-dq 


Then using estimates (4.561 and (4.46) in (4.54) one has 


qt-qst-10 (4.57) 


limsup Ii[p]{t) < C f 1 + Ai(a) + limsup A 2 (t) j . 


Under assumptions A2 the terms limsupA 2 (t) and Ai(a) are bounded. 

t—>-oo 

Thus limsup/i[p](t) < C. Then from formula (4.51) for Ii\p] it follows: 

t—^OO 

(4.58) [ H{\\/p\)dx< [ X(|VTU|)|VfU||Vp|dx+|AQ(t)A.^(t)| tC. 

Ju Ju 

For the first integral in the RHS of (4.58) we have 

(4.59) [ A:(|VfU|)|VfU||Vp|dx < £1 [ H{\\/p\) dx + C\\\7 W\\l.-.^u) + C. 

Ju Ju 

According to Lemma 4.4 in [3] with F 2 {t) as in (4.2) we have 

(4.60) |A^(t)AQ(f)| <£2 [ Hi\Vp\)dx + CF 2 it)\AQ{t)\. 


Choosing £1 and £2 to be small enough and using (4.58) in (4.57) we obtain (4.48). 
Since ^ 2 ( 1 ) is uniformly bounded under assumptions A2, the estimate (4.49) follows. 

□ 


4.4. Estimate of time derivatives. We will obtain the estimate on the difference 
of time derivative of fully transient and PSS solutions. 

Let p(x, t) = p{x, t) — 4>(x, t) and Ps{x, t) = Ps{x, t) — 'i>(ai), then 


eq:trans-shift 


eq:pss-shift 


thr:pt-bound-ah 


(4-61) —= L[p\-'^t{x,t), pIp. =7(i); 

(4.62) ^ = _a = L[WI pjr. =-At. 

Theorem 4.8. Suppose that boundary data satisfies Assumptions AZ-j3. Then for 
any /3 > 2 



Ps,t\^ dx 


|P( + A\d dx < C. 


(4.63) 
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ddt-pt-l-cLh 


Proof. Let Zt =Pt + A and notice that Vzj = Vpj. Subtracting (4.611 and (4.62), 
taking derivative in t, multiplying on \zt\^~^sign{zt), integrating over U and inte¬ 
grating by parts we get 

(4-64) f \ztfdx= [ {L[p]-^SJt-L[W])t\ztf-hignizt)dx 

P dt Ju Jjj 

= [ {V ■{K{\Vp\)Vp))t\ztf-hign{zt)dx- f ^tt\ztf-^signizt)dx 
Ju Ju 

= -(/?-!) [ {K{\Vp\)Vp)fVpt \ztf~'^dx- f ^tt\ztf~^sign(zt) dx 

Ju Ju 

+ j {K{\Vp\)Vp)f N\zt\^~^sign{zt)ds. 

Since -|- ^)|ri = ^j{t) the boundary integral is equal to 
The quantity in the first integral in the RHS of (4.64) is equal to 

-iK{\Vp\)Vp)fVpt\ztf-^ 


= -K{\Vp\)Vpt ■ yptlztf-^ - K'{\Vp\) 


^ (Vpt ■ Vp)(Vp-VpJ ,_ ^_2 


|Vp| 


zt\ 


From (2.13) it follows that 


< aKi\vp\)\vpt\\vptmf-^ 

< aiF(|Vp|)|Vp/|^^|/3-2 + aK{\Wp\)\W'i>t\\Wpt\\ztf-^ 


Then in (4.64) we have 


1 d 
/3 dt 


ddt-pt-la-ah-l (4.65) 


f \ztf dx < - {1 - a){/3 - 1) f K{\Wp\)\Wpt\'^\ztf dx 
Ju Ju 

+ (l + a)(/3-l) [ Ki\Vp\)\V^t\\VPt\\ztf-^dx 

Ju 

+ [ \z,f-^\^u\dx+\iA'^it))^-^Q'{t)\. 

Ju 


To estimate the second term in the RHS of (4.65) we apply Holder inequality and 
Young’s inequality with e = (1-|-a)/(l — a). Since iF(|Vp|) is bounded (see (2.14)) 
applying once more Holder inequality with powers /3/(/3 — 2) and /3/2 we get 


K{\vp\)\v^t\\yptm\^-^ dx < 


1 — a 


Ki\Vp\)\Vp,\^\ztf-^dx 


ineq_l (4.66) 


2(1 + a) Ju 


0-2 

0 


ineq_2 


Similar for the third term in the RHS of (|4.65|) one has 
(4.67) 


J \ztf ^\'i’tt\dx < C\\'Stt\\Lfi{u) (^J \ztfdx + l^. 
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Then in (4.651 we have 


<-C J^K{\Vp\)\Vpt\^\zt\^ dx 


C\mt\\l,^u) ( / %rdx + l 


ddt-pt-la-ali (4.68) 


+ C\\'^tt\\LP(u) (^j dx + + \{A'^{t))d ^Q'{t)\. 


Since the first term in the RHS of ( |4.68 ) is negative, it can be neglected. We 
then have 


z_l 


(4.69) 

where 


^ J \ztfdx<fi{t) J \ztf dx + f 2 {t), 


fi{t) — C'||4'tt||L/3(u) + C'||V4't||^3(y), 

/2(t) = /i(i) + i(A;(t))^-iQ'(t)|. 

The result follows from the Gronwall’s inequality under assumptions A3-/3. 


□ 


In order to prove that limsup Jjj \pf + dx = 0 we will need the Weighted 

i —>-00 _ 

Poincare inequality (see Lemma 2.5., my- 

Lemma 4.9. Let ^ = ^{x) > 0 and function u{x) be defined on U and is vanishing 
on the boundary dU. Assume a > 2. Given two numbers 9 and 9i such that 


> 


(2 - a)* 
there is constant C such that 


and max < 1, 


2n 


710 + 2 


< 01 < 2 - a, 



poincare-w-0 



def:poinc-Mxiu 




thr:pt+A-toO 




ddt-pt-3 


|7t|“dx < C'Mp(C,|u|) 


K{^)\VuY\u\‘^-^ dx 


where 


Mp{^,\u\) = 
92 = 


1+ / dx 

Ju 

0i(0- l)(2-a) 


2 ( 2 -a- 0 i) 


> 0 . 


Theorem 4.10. If Q{t) and 4>(x,t) satisfy assumptions A3-/? and A4, then 

limsup / |p( — Ps jI ^ dx = limsup / \p^ + A\^ dx = t). 
t—^oo Ju ’ t—^oo Ju 

Proof. In (5.15) of [3] we have 

(4.72) 

^ J (Pt + A)'^ dx <-C J A:(|Vp|)|Vptpdx + e 2 J \pt + A\'^ dx + CA3{e2,t), 

where the constant C depends on Jjj iJ(|Vp|) dx and is bounded under assumptions 
A2. 
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pta-0-1 


k—1 


pta-0-2 


We will use Weighted Poincare inequality (4.70) to estimate Jjj K{\\7p\)\\7p^\'^dx 
in terms of Jjj \pf + A\‘^ dx. Namely, let u =%+ A — and ^(x) = |Vp|. Then 

mIt; = 0 and Vit = Then we have 

(4.73) [ \pt+A\‘^dx= [ \u + dx < ^ [ \u\'^dx+^ [ \A'^{t)\'^ dx 

Ju Ju ^ Ju ^ Ju 


< CMp{^, |u|) 


K{Wp)\Wp,\^dx 


+ ^|A;(t)p. 


Here 


Mpi\Vp\,\u\) = 


1+ / |Vpp-“ + |p*+H-A;(t)r dT 


< c 


in view of Theorem |4.7| under Assumptions A2 and Theorem |4.8| under Assumptions 

A3-13. 

Thus 


|P( + Ap dx < C 


K{Wp)\Wpt\^ dx 


\U\ 


|A;(t)r 


On the other hand note that qtdx = 0 and gdr, = {qt — i?s,t)|ri = A!(,(t) + F{{t). 
Here F[{t) = Aq)^) + J^jA^tdx (see (|4.1|)). Notice that limsup |F((<)| = 0 under 

t—>-oo 

conditions A4. Then by Trace theorem, Poincare inequality and the fact that 
Jjj qtdx = 0 we have 

|A;(t)p<c(^^ |g*|ds + C'|F((t)|^ <c(^J^\Vqt\dx'^ +C\F[it)\^. 

Next, as Vg* = Vpj, applying Holder inequality we have 

(4.74) (^J^\Vqt\dx'^ <C J^Ki\Vp\)\Vpt\^dx-J^K{\Vp\)-^dx. 

Notice that from (2.14) and Holder inequality it follows 

(4.75) [ K{\\7p\)-Ux < C [ (1 +|Vp|)“da; < C + C [ |Vpp-“ dx < C, 

Ju Ju Ju 

The last estimate in (|4.75 ) follows from Theorem 4.7 under conditions A2. 

Thus 

|A;(t)|2<c/ K{\Vp\)m\^dx + C\F{it)\^. 

Ju 


Then in (4.73) we have 


(4.76) ^ \pt + Ap dx<C (j® {t) + m) + C\F[{t)\^ = C ■ p(l{t) + \F[{t )\^), 


where 


(4.77) (^(s) = sHs + |F((t)|2; I{t)= [ Ki\Vp\)\Vptfdx. 

Ju 

Due to condition on |A((t)| 



















WELL PRODUCTIVITY INDEX FOR COMPRESSIBLE FLUIDS AND GASES 


22 


ddt-pt-4-new 


th:pi-pipss-0 


sec:dirichlet 


dir:dl 


Then due to definition (4.6) of ^ 3 (^ 2 ,^) one has in (4.72) 

(4.78) 

^ J + A)^ dx < - Cip~^ \p^ + A\'^ dx'^ j dx + CA3{e2,t)- 


Let y{t) = fjj \pf + A\^ dx. Due to Theorem ( |4.8| ) one has in (4.78) 

(4.79) y'(t) < -Cfivit)) + CA^{e2,t) + ^£ 2 , 
where 

(4.80) /(s) = y)“^(s), /(s) > 0 for s > 0, and /(O) = 0. 


Under the Assumptions A3-/3 we have 0 < y(t) < C. According to Lemma 4.5 


we 


will get limsupt_^oo y{t) < C'(/j(limsup(_^oo(^ 3 (^ 2 ,i)+C 2 )- Under Assumptions A4 
limsupj_^oo A 3 (£ 2,0 = 0 ^od therefore 

limsupj/(t) < p{e 2 )- 


Taking in above £2 —>■ 0 we complete proof of the Theorem |4.10[ 


□ 


Finally Theorems 4.6 4.7 and 


A1-A5 allow us to conclude that 


4.10 and Lemma 4.1, [3], under the conditions 
|V(p-Ps)|| 1 ^ 2 -a —>■ 0 as t —> 00 in Lemma |4.1 


Combining it with Lemma we get the main result for the solutions of IB VP-1: 
Theorem 4.11. If boundary data Q{t) and 4>(a;,t) satisfy assumptions A5 then 
(4.81) Jgit) - Jg,PSS 0 as t —>■ 00 . 


5. lBVP-11: Asymptotic convergence of Diffusive Capacity 


Let p{x,t) be the solution of IBVP-Il (3.1), ( |3.3| ), (3.4) with Dirichlet boundary 
condition ( |3.8[ ) with boundary data 'l'(a:, t). Let Jg{t) be the corresponding diffusive 
capacity. Let ps{x,t) = — At -I- lU(a:), lU(x)|ri = </ 2 (x), be the PSS solution of this 
IBVP defined by the boundary data —At, A = Qs/\U\, and $(a;). Let Jg,pss be 
the corresponding diffusive capacity. 

5.1. Assumptions on the boundary data. Let 

(5.1) D{l3,t) = \\'^tt\\LP{u) + \\'^'^t\\lp(^u) for any/? >2. 

We also will use some notations from [13]. For a > b = 


To = 


"(2-a) Ipt 

(2-a)(ra+l)-n’ 


(5.2) 


(5.3) 


A{a,t) = \\Vnir 


G4{t)= / 

Ju 




Vl-a) 


ro(l-a) 


JU 
^ dx 


14 - 


tllL'-o 


\^tt\ dx. 


+ / |V4't|^dx-k / \'Ut\'^dx + 

Ju Ju Ju 

The results in this section are obtained under the following conditions on the 
boundary data 
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assumpDl 


assumpD2-l 


assumpD2 


assumpD3-l 


assumpD3-2 


eq; dir-trans-shift 


eq:dir-pss-shift 


thrm:pt-ah-bound-dir 


ddt-pt-dir 


/ C30 

£)(/3, t)dt < oo for any (3 >2. 


Assumptions Dl: 

(5.4) 

Assumptions D2: 

(5.5) limsup (A(a, t) + G 4 (t)) < C; 

t—>-oo 

(5.6) limsup ( Z?(2, t) + j + Ap da; | = 0. 

t^OO \ JlJ J 

Assumptions D3: 


(5.7) 

(5.8) 


limsup A(a,t) + II < C'; 

t—>oo 

limsup||V(A4,)||i2([;) = 0. 


According to Lemma 3.2 in order to prove Jg{t) —>■ Jg,pss as t —>■ 0 we need to 
prove the following estimates: 

(5.9) ^q(^) 0 as t ^ oo see Lemma [5?2{ 

(5.10) / |V(p — Ps)p““ da; — > 0 as t ^ oo seeTheorem 15.31 
Ju 

5.2. Estimates on the solution of IBVP-II. Along with the solutions p and Ps 
we will use the shifts p{x, t) = p{x, t) — 4'(a;, t) and Ps{x, t) = Ps{x, t) — $(x) + At. 
These shifts Obvioulsly satisfy 

dp 


— = L[p] - 4't(a;,t), p|r = 0; 


(5.11) 

(5.12) 

Lemma 5.1. Suppose boundary data '^{x,t) satisfies assumptions Dl, then 


^=0 = L[W]+A, p,|r = 0. 


(5.13) 


/ dx < C for any fi >2. 

Ju 


Proof. Similar to (4.64) we have 


(5.14) 


< - (1 - a)(/3-1) J K{\Wp\)\Wpt\'^\ptf ^ da; 

+ (l + a)(/3-l) [ Ki\Vp\)\V^t\\^pt\\pt\^-^dx 
JU 

+ [ \ptf-^\^tt\dx. 

Ju 


Similar to (|4.66|) - (|4.69|) we have 

|pt|^ da; < Il(^,t) / \ptf dx + D{j3,t). 


(5.15) 


d 

dt 


The result follows from the Gronwall’s inequality under assumptions Dl. 


□ 
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dir: int_li 


lem:deltaq-0 


est:pta2 


delta-q2 


pbart-a-1 


pbart-a-2 


pta5 


Further we will need the following result obtained in |13] : under assumptions 
D2.1 ra (see [B], Th. 4.5) 


(5.16) 


[ |Vpp-“da; < C / H{\Vp\)dx + C <C. 

Ju Ju 


Lemma 5.2. Suppose boundary data satisfies assumptions D1 and D2. 

Then 

(5.17) limsup MAQ(t)| + (7 / |pt + Ap da; | = 0 as t —>■ oo. 

(->•00 \ Ju ) 

Proof. First, notice that 

d f d f 

— J p{x,t) dx = Q{t) and — J ps(x,t) dx = Qs- 


Thus 


\AQ{t)\^ = \Q{t)-Qs\^ = 


Pt dx - Qt 


/ {pt + A) da 

Ju 


(5.18) — ^ / \pt+A\^dx<C f Iptl"^ dx + C f \'^t+A\'^dx. 

Ju Ju Ju 

Assuming + Ap da; —)■ 0 as t —)■ oo it is sufficient to prove that dx ^ 0 

as t —)■ oo. 

Differntiating both sides of (5.11) in t, multiplying by pt and integrating over U 
we have due to zero boundary condition 

(5.19) J p^dx = -J (Ar(|Vp|)Vp)tVpt da; - y ^ttPtdx. 

Similar to estimate (5.10) in [3] we get since Vp = Vp — 

(5.20) 4 f ^tdx<-C [ K{\Wp\)\\/pt\‘^dx + e( [ \pt\‘^dxy~^ +CD{2,t). 


dt 


where 

(5.21) 


Di2,t) = \\^tt\\ir 




dlL2- 


Applying the Weighted Poincare inequality (4.70) with u = pt, ptlui =0 and 
^ = Vp with powers a = 2 and 9 = the first integral on the RHS of (5.20) 

we get: 


-I Pidx<-CMf; 


-9 


IPtpda; 


+ £ 


IPtpda; 


CD{2,t). 


Here 

(5.22) 


Mpi\Vp\,\pt\) = 


1 + / |Vp|^-“ + \pt 


|26»2 


dx 


< C 


by virtue of (5.16) under assumptions D2.1 and Lemma 5.1 under assumptions Dl. 
Choosing sufficiently small e we get 


(5.23) 


limsup / |ptp da; < limsup Zi(2, t). 

t^oo Ju t^OO 


Under assumptions D2.2 the RHS of (5.23) converges to 0 as t —oo and the result 
of Lemma follows from (5.18). □ 
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dir: pbax-psbcir 


th:dir:p-ps 


phi-2-a 


th:pi-pipss-dir 


Further we will use another result from m- under assumptions D3.1 ( |5.7[ ) (see 
[13], Th. 4.3) 


(5.24) 


[ \p{x,t)\'^+ \psix, t)\^dx <C. 
Ju 


Theorem 5.3. Suppose the boundary data 4'(x,t) and $(x) satisfies assumptions 
D1 — D3. Then 


(5.25) 


l|V(p-p.)|| ^ 2 ——y 0 as t —y oo. 


Proof. From (2.16) we get 
(5.26) 


\Vip-ps)r‘^dx] <C^ / {Ki\Vp\)Vp-Ki\Vps\)Vps)-V{p-ps)dx, 


where (7$ = C,j,{p,ps) is as in (2.17). Since V(p —Ps) = V(p— W) + V(A,p), (7$ 
is bounded by virtue of (5.16) and (4.23) under condition ||V(A^)|ji 2 -a([/) < C. 
Since p — ps\ri = 0, integration by parts in the RHS of (5.26) gives 

W^ip- Ps)\\l2-a <- C [ (pt- Ps,t)ip-P 8 )dx 

Ju 

+ C [ iKi\Vp\)Vp-Ki\Vps\)Vps)-V{A^)dx. 

Ju 

Since Kf) is bounded we have 


(5.27) 


(i7(|Vp|)Vp - K{\Wps\)Wps)WiA^)dx 

Ki\S/p\)\Wp\^dxY ( [ Ki\Wp\)\ViA^)\^dx 


< 


Ki\VW\)\VW\^ dx 


Ki\VW\)\ViA^)\‘^ dx 


< 


Hi\Wp\)dx] + / |VVF|2-“dx 


|V(A,f)|ji2([7) 


< C||V(A,f)|ji2([7). 


The last inequality holds in view of (5.16) and (4.23). 

Thus 

W'^(p - Ps)\\l2-a < - C [ ipt +A)ip-ps)dx + C\\\'iA^)\\L^u}■ 

Ju 

The result follows in view of Lemma 5.2 and estimate (5.24|) under the assumptions 
D1 - D3. 

□ 


Lemma |5.2| and Theorem |5.3| combined with Lemma |3.2| prove the main result 
for the solution of IBVP-11: 

Theorem 5.4. Suppose the boundary data satisfies assumptions D1 — D3. Then 
Jgit) — Jg, PSS 0 as t —>■ oo. 
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sec:Gas 


6. Productivity index concept for the flow of ideal gas. 

In this section we discuss the concept of the Diffusive Capacity for an ideal gas 
flow in the porous media for a well-reservoir system. We will repeat some discussion 
from Sec. and provide more details. The equation of state for an ideal gas takes 
the form (see [10], [18]) 


eq: state-gas (6.1) 


p{p) = Mp, 


where p is the density of the fluid and M is a proportionality constant. Without 
loss of generality we assume M = 1. As a momentum equation we consider the 
second order Forchheimer equation (2.1) which takes the form (see, for example, 

m) 


eq:2forch-l 


( 6 . 2 ) 


au + Ppu\u\ =—Vp, j3 = F^k 


where F is the Forchheimer coefficient, $ is porosity and k is permeability of the 
porous media. 

As in case of slightly compressible fluid, the above equation can be solved for u 
and rewritten in terms of a nonlinear permeability function K 2 {\'^p\)' 


eq:2forcli (6.3) 


u = —K 2 {p\Vp\)Vp = — 


a -I- + 4/?p|Vp| 


Vp. 


Combining (6.1) and (6.3) together with the continuity equation (2.6) yields the 
parabolic equation for pressure 


gas-equation 


(6.4) 


dp 


L[p] = V • (iF2(p|Vp|)pVp) = V • 


2p 


a + + 4/3p|Vp| 



As before the domain U models the reservoir with boundary split in two parts Fj 
and Fe. F^ is the well-boundary while Fg is an exterior non-permeable boundary of 
the reservoir. 

We impose Dirichlet Pseudo Steady State (PSS) boundary conditions on F^ 


bc-gas:dirichlet 



bc-gas:neum 



gas-init-cond 


(6.5) p|r. =B-At 
and zero mass flux boundary conditions on Fg 

(6.6) pu ■ A^lr, = 0 
and the initial condition 

(6.7) pix,0) = 1/52 +(j)o{x) 
for (l)o{x) > 0. 


Remark 6.1. Constant R is a positive parameter (generally large) characterizing 
the initial reserves of gas in the reservoir domain U. The constant A is associated 
with the amount of gas extracted at the well-bore F^. Thus the quantity B — At 
quantifies the gas reserves in the reservoir at the moment t. 


For gas filtration in porous media, engineers define the Productivity Index as, 
see [51|7]: 
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Pl-gas-def 


gas-total-flux (6.8) 


def:pi-gas (6.9) 


eq:gas-W (6.11) 


eq:gas-Wb 


eq:gas-Wc 


f luxW 


gradW 


gas:eq:pO 


gas:eq:pOc 


Definition 6.2. Letp{x,t) be a classical solution of equation 6.4 satisfying bound¬ 
ary conditions 6.5 and \6.(^ Let Q(t) he the the total mass flow through the well-bore 
boundary 


Qlf) = J pu ■ Nds = J K 2 {p\Vp\)pVp ■ Nds, 


where N is an outward normal toTi. Then productivity index for the well-reservoir 
system is defined by 


J = 


Q{t) 




As we saw earlier, in case of slightly compressible flow under the boundary 


conditions (6.5) and (6.6) the value of the PI, defined as in (2.10) is stabilizing to a 


constant value (3.14), which is defined by the PSS solution and depends only on the 


paramerter A. Obvioulsy, in case of the gas flow descirbed by (6.4) this feature is 


not valid anymore and the PI is not stabilizing in time. Nevertheless the approach 
developed in previous sections is applicable to the gas flow as well. Namely we will 
introduce special auxiliary pressure function characterized by the time independent 
PI. We will then study the relation between the general time dependent PI and this 
time independent one. 

Let Qo = be a given constant mass-rate of gas production and consider the 
auxiliary pressure distribution given by 


p_0def (6.10) 


Po 


(x,t) = yJ{B - Aty +2W{x), 


where W {x) is the solution of the following BVP 


- V • {K2{\VW\)VW) = -V 


2VW 


a + +AI3\VW\ 


= A, 


( 6 . 12 ) 

(6.13) 


W = 0, on 
VW-N = 0, 


r*, 

on 


Integrating Eq. (6.11) over the domain U, using integration by parts and boundary 


condition (6.13) yields to the integral flux condition 


(6.14) 


J -K2i\\7W\)VW-Nds=\U\A = Qo. 


From (6.10) we can also derive that 
(6.15) 


PoVpo = VVF. 


Combining this property, definition (6.10) and BVP (6.11 )-(6.13) it can be easily 


verified that the auxiliary pressure po satisfies the following equations 

dpo 


(6.16) 

(6.17) 

(6.18) 


dt 


- L[Po] = foix,t) in C/, 


Po{xA) = B — At on Pi, 
PoVpo ■ = 0 on Pg. 
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gas:def-fO 


fluxpO 


jpo 


time-constr-1 


time-constr-2 


where 

(6.19) 


foix,t) = ^ - V • (i^ 2 (po|Vpo|)poVpo) 


-A{B - At) 

+W 

-A{B - At) 


- V • iK2i\yW\)VW) 
+ A 


y/{B - At)^ + 2W 

- 4 fi B-At \ ^ ^ 

“ y ^/{B - At)^ + 2 w) ^ 


From (6.14) and (|6.15 1, it also follows the total mass flux condition for po 


( 6 . 20 ) 


^ -K 2 ipo\Vpo\)poypo ■ Nds = -K 2 {\VW\)VW ■ Nds = Qo. 


Using this last result and the explicit representation of po{x,t) in formula (6.9) 
leads to the following Proposition. 

Proposition 6.3. The Productivity Index for gas flow defined on po{x,t) is time 
independent and is given by 


( 6 . 21 ) 


J[Po\ = 


Qo 


j^fu2Wix)dx' 


Numerical computations, performed for several basic reservoir geometries, show 


that if the initial data in the system (|6.3[)-(6.6 ) is given by p(a;,0) = po(a^,0), then 


the corresponding productivity indices J[po] and J[p\(t) are almost identical for a 
long time, see Fig. as long as the quantity 

, , B- At 

( 6 . 22 ) 


V'(B - At)"^ + 2W 


1 , 


or, equivalently, as long as 
(6.23) {B-Atf:^2\\W\\c 


t ^ Tcrit ^ 




A 


In case of compressible flow, we have to fix a critical time 

T 

crit — ^ I 


and for t > Tcrit the negative boundary data (6.5) leads to the violation of ellipticity. 
This behavior is qualitatively justified by the fact that as long as ( |6.22[ ) holds, 
function fo{x,t) in the RHS of ( |6.16 l is negligible and p{x,t) and po{x,t) behave 
similarly. On the other hand, as (B — At) —)■ 0 then fo{x, t) approaches the constant 
value A, and the two solutions diverge from each other. Then two productivity 
indices J[po] and J[p]{t) also diverge from each other. 

This phenomenon is observed on the actual field data and has a clear practical 


explanation. Notice that the denominator in the formula for J[po] (6.21) 




\U\ 
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Figure 1. Comparison between the time dependent productivity 
index J\p\ and the PSS productivity index J[po] as the time t 
approaches the critical value T^rit = ^ = 2000. 


fig:PIgas 


(since W > 0 for any x S C/) is a measure of the pressure drawdown needed to 
maintain constant production Qq. As long as the gas reserves are considerably 
larger than the the pressure drawdown 

(6.24) ^B-Atf>C>2\\W\\^>^^^, 

then the distributed source term /o (equivalent to reservoir fluid injection) needed 
to maintain constant production Qq is negligible. Otherwise when the gas reserves 
are comparable in magnitude with the pressure drawdown, then a possible way to 
maintain constant production rate is by resupplying the reservoir by fluid injections. 

In the remaining part of this section we will theoretically investigate the dif¬ 
ference between the functions p{x, t) (the actual solution of the problem) and 
Po^Xjt) (the solution of the auxiliary problem) depending on the key parameter 
Tcrit = B/A. Unfortunately for the Forchheimer case we were not yet able to ob¬ 
tain the appropriate estimates for the differences between pQ{x,t) and p{x,t). We 
will report mathematically rigorous result only for the case of compressible Darcy 
flow. We will show that for the fixed Tq for all times t G [0, Tq] when p{x, t) > 0 
on U X [0,To] the productivity indices J[p] and J[po] are becoming closer to each 
other with the increasing parameter B. The obtained results make us believe this 
comparison can be proved in the general Forchheimer case as well. 

We consider the case of positive solutions. For that assume that our time space 
domain belongs to the time layer 0 < t < Tq < Tcrit- Let D = U x (0,To]. The 
following useful inequality follows directly from maximum principle 


Lemma 6.4. Let p{x, <) > 0 for all 0 < t < Tq < Tcrit is a classical solution of the 
equation (6.4) with boundary conditions (6.5) and (|6.6[). Then for any 0 < t <Tq 


(6.25) 


p>B-At 


p{x, t) > min{p(x, 0); i? — At}. 
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lem:comparison 


L-operator 


u-bar-eq 


u-eq 


neumann-on-gaiiraia2 


bc-comparison 


comparison 


lemma-identity 


gas :p-pO-identity 


Proof. Inequality (6.25) follows from the maximum principle (see, for example, 
m), since the nonlinear equation ( |6.4[ ) is uniformly parabolic in the domain D, 
and the boundary function in (6.5) is decreasing with time. □ 


6.1. Analytical comparison of the solutions for case of the Darcy flow. 

First the following maximum principle follows from the results in |22j . 


Lemma 6.5. Let (x,t) C D = U x (0,T], U C 
defined on U x (0,T] 


(6.26) 


/I ^ ^ ^i,j (^j 






dxjdxn 


'^bi{x,t) 


and an elliptic operator C is 

d 


dxi 


with C 1^1 Ei=i \bi\ < C. 

Let dU = Fi U F 2 , where for simplicity Fi and r 2 are nonintersecting compact 
sets. Let dD = U x {0} U dlJ x (0,r] be the parabolic boundary of the domain. 

Assume u{x,t) and u{x,t) > 0 in D are the solution of inequality (6.27) and 
equation (|6j2|) correspondingly: 


(6.27) 

(6.28) 




Assume also that u(x, t) and u(x, t) satisfy the homogeneous Neumman conditions 
on the boundary F 2 : 


(6.29) 


du 

dn 


r2 


du 

dn 


= 0 . 


r2 


The following comparison principle holds: if 

(6.30) u> u on U X {0} U Fi x (0, T], 
then 

(6.31) u > u in D. 

We will now obtain some integral comparison results between actual and auxil¬ 
iary pressures in case of linear Darcy flow. Namely, let the Forchheimer coefficient 


F = 0 in (6.2). Let po{x^t) be the auxiliary pressure given in Eq. (6.10) andp{x,t) 


be the classical solution of IBVP (6.4) with the PSS boundary conditions (6.5)-(6.6) 


and initial data given by p{x, 0) = Po{x, 0) = P 2W{x). First, we will prove a 
useful integral identity for the difference p{x, t) — po{x, t). 

Lemma 6.6. Suppose T < T^rit = BjA. Then the following identity holds 

pT p ^ r / \ 

(6.32) 


/7 

/o Ju 


iP + Po)ip-pofdxdt+^J^J^ {p‘^-pl)d?j dx 
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where fo is defined in (6.191. The identity above can be further rewritten as 


gas:p-pO-identity-l (6.33) 


f\ 

/o Ju 


iP + Po)ip-Pofdxdt+^J^J^ {p^-pl)d?j dx 

I I / fo{x,T)dT^ dtdx . 


Proof. For the Darcy case the equations for p{x,t) and pQ(x,t) take the form 


gas:darcy-p 


(6.34) 


gas:darcy-pO 


(6.35) 


^ = V • (pVp) = ^A(p2), 

^ = V • (poVpo) + foix,t) = ^A{pI) + fo{x,t). 


Subtracting the second equation from the first one yields 

^(P-Po) = _ fo{xfi). 

Following the ideas of Oleinik (see [H]), we multiply both sides by the test function 
{p^ — Po) dr, integrate in time form 0 to T and in space over the domain U. 
Then it follows 


eq:miilt-by-test (6.36) 


dtdx 

- j J ^/o ^ (P^ - Po) dr'^ dt dx. 


By using integration by part and the fact that p{x, 0)—po(a:) 0) = 0 and {■)dt = 0, 
the integral on the LHS of (6.36) can be rewritten as 


est:gas-lhs 


(6.37) J J ^^(p-po)^ {p'^-pl)dT^ dtdx 

= -J J ~ {p^ - pI) d'^ dt dx 


{p-Po){p'^-Pl)dtdx= [ [ {p + po){p-Pof dtdx. 
u Jo Ju Jo 
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By using the divergence theorem the first integral in the RHS of (6.36) can be 
rewritten as 


est:gas-rhs (6.38) 


^lulo 

= -ll£ (v(/-pg)-/"v(p^-pg) 

^ '^iP^~Po)dT\ dtdx-^J^ V{p^-pl)dT 


— pi) dr dt dx 


dx. 


Substituting (6.37) and (6.38) back in (6.36) we obtain (6.32). Finally in order to 
obtain the alternative identity (6.33) we use the following integration by part for 
the right hand side of Eq. (6.32) 

lu Jo 


I I {v^{^^'^)-pli.^^'^))dTjdtdx 

J fo{x,T)dT {p‘^{x,t) - pl{x,T))dT 


T / pt 


/O \^0 
pT 


0 / 
dx 


/o(a;,r)drj [-{p‘^{x,t) - pl{x,t))] dt 

lu I Mx,T)dT^ dtdx . 


□ 


From the above Lemma follows 

Proposition 6.7. Under the conditions of Lemma \ 6. d| then the following compar¬ 
ison holds 


gas:p-p02-f02 (6.39) 

T 


gas:p-pO-estimate (6.40) 


[ [ {p^—pi)dxdT < C / f [p^—pifjdr dx < C f [ f fo{x,t)dt\ dx. 
Jo Ju Ju Jo Ju \Jo ) 

Proof. Since fo{x,t) > 0 for all {x,t) (see (6.19)) and according to Lemma 6.5 
pKxjt) — P^ix.t) > 0 for all {x,t), in the RHS of (6.33) we have 

'X' / i \ 'X' 7^ 

/ f \{Po-P^) f fo{x,T)dT] dtdx < [ [ (Po-P^)/ fo{x,T)dTdtdx 

Ju Jo \ Jo J Ju Jo Jo 

= f f {Po-p'^)d* [ fo{x,t)dtdx 
JuJo Jo 

<£ J (pI-P^) dt^ dx + C^ J fo{x, t ), 


I dt dx. 
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gas-poincare 


constr-deriv 


constr-pt 


gas:p-pOtoO 


ueq 


uBC 


a-estimate 


By Poincare inequality we have 
(6.41) 


-I V I- -I V I- 

f ( - Po) <C [ [ -pI) dr dx<C [ V / -pl)dT 

Jo Ju Ju Jo Ju Jo 


dx. 


Finally, estimating the RHS of (6.33) using (6.40), neglecting the first term in 
LHS, and using (6.41) for the second one, we get: 

2 


T ^ / T f T 

[ [ {p^-Po)dx dr <e [ i [ {pi - p'^) dt\ dx + [ ( [ fo{x,t) 

Jo Ju \ Ju \Jo J Ju \Jo 


Choosing appropriate e we get (6.39). 


dtj dx. 

□ 


6.2. Stability of PI with respect to initial and boudnary data. In this 
section we will show that for the given time interval 0 < t < Tg < T^rit the difference 
between the Pi’s for the actual and the auxiliary problems becomes small as the 
parameter B for the initial reserves becomes large. As before we consider only the 
linear Darcy case with F = 0 in (6.2). 


Let p{x,t) be the classical solution of IBVP (6.4)-(6.6) with the corresponding 
productivity index J[p]{t) defined by (6.9). Let pQ{x,t) be the auxiliary pressure 
given by (6.10) with the corresponding productivity index J[po\ defined by (6.21). 
In linear case the original pressure p{x,t) inherits the following properties of the 
auxiliary pressure po{x,t) (for the details see [12]): for alH, 0 < t < Tg < Tcrit 

(6.42) \Vp{x,t)\ + Ap{x,t) < C < oo] 

(6.43) \pt{x,t)\ < C < oo. 

In the following two lemmas we will obtain estimates for p — po and pt — po,t as 
B —>■ 0. These will be used later to prove Theorem |6.10| 


Lemma 6.8. Under the constraints (6.42) for B > 1 

C 


(6.44) 

Moreover, if B 

(6.45) 


0>P-PO>- 




and 


0>p^-pI>- 


p — pg —>■ 0 and p^ — 


pI 


0 . 


Proof. Let a{x,t) = p + pg, bi{x,t) = d{p + po)/dxi, i = l,...,n and c{x,t) = 
A(p + pc). The direct calculations show that the function m = p — pg is a solution 
of 


(6.46) 

(6.47) 


Ut — aAu + biUx, + cu = —fo{x, t) in U x (0, Tg], 


i=l 


= 0 on dU X (0,Tg]. 


Under the conditions (6.42) it follows that \bi{x,t)\ + c{x,t) < C independently of 
time t. Eq. (6.46) is parabolic in it = p — pg and 


(6.48) 


Cii? > a{x,t) > CqB 
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for some Co,Ci > 0. Thus following the standard arguments using the barrier 
functions (see m) we get that 


estimate p-pO (6.49) 


— C max foix^t) < p — po < 0. 

<To 


Taking Tq < \Tcrit = one can get that 


/o< 


A 


< 


C 


VB^ + 2W+^B2 - 1 + B2 


0 as i? —)■ oo. 


Then the result follows from (6.49) and (6.48). 


□ 


gas:pt-pOt-toO 

Lemma 6.9. Under the constraints ( 

6.42 

) and ( 

6.43 


(6.50) 


Pt — Po t ^ 0 HS B ^ oo. 


Proof. Multiplying equations (6.34) and (6.35) on p and po correspondingly we get 

2ppt-pA{p'^) =0, 

2poPo.t - PoA{pI) = fopo . 

Let u = ppt = ^ip^)t and uq = poPo,t = ^(Po)t- Then differentiating the equations 
above in t we get 

2ut -pAu = ptA{p ^), 

2wo,t - PoAuq = po.tA( pq) + 2/o,tPo + 2/oPo.t • 

Since 2pt = A{p^) and 2po.i — /o = A(po) we can rewrite 
2ut — pAu = 2p( , 

2Mo,t - PoAuq = 2pl t + 2/o,tPo + foPo,t ■ 

Subtracting two last equations from each other we get 


eq:u-u0 (6.51) 


d 

‘^dt "“o) - pAu + poAuo = 2{pl-pl ^) - Fi{x,t) 


where Fi(x, t) = 2fQ^tPo + foPo.t- Denoting z = u — uq and adding on both sides of 
(6.51) the term pAuQ we get 

dz 

2— -pAz = 2 {pI-pI ^) - {po-p)Auo - Fi{x,t)- 
Further, since pt = u/p and po,t = uo/po we have 


eq:z-l (6.52) 


^dz 

^dt 


U Uq 


P Po 


-pAz = 2- ] {pt +Po,t) - iPo-p)Auo - Fi{x,t)- 


In the first term in RHS adding and subtracting the term uqPo in the numerator 
we get 

U Uq _ Pou - UoP _ Poiu - Uo) Uo(j) - Po) _ Z {p - Po) ^ 

P 


Po 


PPo 


PPo 


PPo 


P 


P 


-Po,t- 


Thus from (6.52) we have linear equation for z 
dz 


2^ -pAz = 2-C{x,t) +F 2 {x,t) - Fi{x,t), 

ot p 
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where 


C{x,t) = {pt+Po,t) 


2po,t 


F 2 {x,t) = {p-po) Amo- -{pt +Po,t) 


P 


Under conditions (6.42) and (6.43) and in view of Lemma [6^ \C{x,t)\ < C < oo 
and F 2 {x,t) +Fi{x,t) —)• 0 as i? —)■ oo. Then \z\ 0 as B 

since 


00. The result follows 
z = u-uo=ppt- PoPo,t = Po{pt - Po,t) - Pt{Po - p)- 


□ 


Finally, we will state the main theorem of this section. 



theo:gas-pi 

Theorem 6.10. Under constraints (6.42 

) and ( 

6.43 

) 

gas:cond-compar-pi 

(6.53) J[p]it) ~ J[Po\ “t 0 for B ^ oo 


Proof. From the boundary condition plr^ = B — At and the fact that {B — At)'^ = 
pl{x,t) — ‘2W{x) (see (6.10)) we get 

1—7) ](A — _ Qi^) _ I _Qo_ 

— /p pVp ■ N ds ~ fr Po^Po ■ ds 
jF p2 dx - (B - At)^ jFf^2W(x)dx 


ifuP{x,t) dx 


dt Ju 


1^ lu^P^ - pI) dx + Jir\ (^) 


Adding and subtracting the term po{x,t) dx in the numerator of the first 
fraction we get 

{j[po]-j[pm)-\u\ = 

d 


TtIu(P-Po)dx 


jFfu2W(x) dx 
le numerator i 

Ju fo(,x,t) dx 


+ ^| / 


luiP^ — Pq) dx + 2 J^W dx 2 JjjW dx 

1 1 


fjj(p^ — Pq) dx + 2 JjjW dx 2 JjjW dx 


Tt!u^P-Po)dx 


i iluPodx) 


- luiP"^ - Po) 


Ju(p^ — Pq) dx + 2 J^W dx 2 JjjW dx JjjiP^ ~ Pq) dx + 2 f^Wdx 

fufojxA) dx 


2 J^Wdx 


Finally 

{j[po]-j[pm)-\u\ = 


Jufo{x,t)dx 


+ 


jp-jpO (6.54) 


2 JjjW dx Jjj {p^ — Pq) dx + 2 f^W dx 
1 1 (p - p„) dp - / ^ - pD dp ■ 


Since since /q —>■ 0 as B —>■ oo (see (6.19)) then 

Jufo{x,t)dx 


2 JjjW (x) dx 


0 as B ^ oo. 
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Then (6.53) follows from Lemma 6.8 and Lemma 6.9 


□ 


7. Conclusions 

• The notion of diffusive capacity/productivity index is studied in case of 
Forchheimer flow of slightly compressible and strongly compressible fluids. 
In general case the PI is a time dependent integral functional over the 
pressure function. 

• In case of slightly compressible fluid we consider two types of boundary pro- 

hle on the well-boundary: the total flux condition and Dirichlet boundary 
condition. In both cases we prove the convergence of the time dependent PI 
to a constant value without any constraints on the degree of nonlinearity of 
Forchheimer polynomial. This generalizes our previous work [3], requiring 
the estimates on the mixed term a yf 2, and resulting in 

stronger constraints on smoothness of boundary data. 

• In case of strongly compressible fluid, the ideal gas, we study the Dirichlet 
boundary problem. The quantity B — At prescribed on the well-boundary 
specifies the gas reserves at the moment t. We show numerically that the 
PI stays the same until it suddenly blows up when time approaches the 
critical value T^rit = B/A. This fact corresponds with field observations. 
We associate the constant PI with the special pressure po{x,t) and in case 
of linear flow of the ideal gas we analytically study the relation between 
general pressure and Pq{x, t). The results on stability of the PI with respect 
to the initial gas reserves are obtained. 
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